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We compute the phenomenological signatures of a model (Watanabe et al' 09) of anisotropic 
inflation driven by a scalar and a vector field. The action for the vector is U(l) invariant, and the 
model is free of ghost instabilities. A suitable coupling of the scalar to the kinetic term of the vector 
allows for a slow roll evolution of the vector vev, and hence for a prolonged anisotropic expansion; 
this provides a counter example to the cosmic no hair conjecture. We compute the nonvanishing 
two point correlation functions between physical modes of the system, and express them in terms of 
power spectra with angular dependence. The anisotropy parameter g, for the scalar-scalar spectrum 
(defined as in the Ackerman et al '07 parametrization) turns out to be negative in the simplest 
^s^j . realization of the model, which, therefore, cannot account for the angular dependence emerged in 

some analyses of the WMAP data. A g* of order -0.1 is achieved when the energy of the vector is 
about 6—7 orders of magnitude smaller than that of the scalar during inflation. For such values of the 
parameters, the scalar-tensor correlation (which is in principle a distinctive signature of anisotropic 
spaces) is smaller than the tensor-tensor correlation. 



in 

i ,' I. INTRODUCTION 

o. 

The WMAP data on the CMB anisotropies [1J have strongly improved our knowledge of the Universe, and are in 
overall excellent agreement with the theory of inflation. However, the unprecedented quality of the data has driven a 
number of analysis on finer effects, which seem to be hard to reconcile with the simplest inflationary models. These so 



hemispheres [7|], and broken rotational invariance [8]. Most of these effects (and, particularly, those appearing at 
the largest scales) can be generally thought as violation of statistical isotropy of the CMB. It has been suggested in 
@ that a period of anisotropic expansion during inflation may explain some of these features. Specifically, ref. @ 
pointed out that the anisotropic expansion would provide a primordial power spectrum with a dependence on the 
directionality of the modes, and computed the corresponding correlation matrix (ai m a^, m ,) tfcSu' 5 mm i. 
Ref. 10] introduced the power spectrum parametrization 

P(k)=P(fc)[l + S.(*K 2 ] (!) 

where £ is the cosine of the angle between k and a given fixed direction, and k = |k|, and computed the corresponding 
CMB correlation matrix as an expansion series in g*. Althou gh the specific inflationary model proposed in (Toj to 
obtain this power spectrum was later found to be unstable [iTl.ll2j , the ACW study based on (JTJ is extremely useful, 
since it provides a general reference result for generic breaking of rotational invariance during inflation. 1 The ACW 
parametrization was tested against the WMAP data in Hi, where a 3.8<r evidence was obtained for nonvanishing 
(an upper limit on the anisotropy was instead obtained in |15j|). The study assumed a constant g*, and can be likely 



be applied to slow roll inflationary backgrounds, in which a mild scale dependence may be expected. This analysis 
was then refined in two more recent works [l6l.[l~7j. which include a (— i) 1 factor in the covariance matrix which was 
neglected in the first version of [loj ]. and in the analysis of Q. The inclusion of this factor has two important effects. 
Firstly, it substantially increased the significance of a nonvanishing asymmetry: g» = 0.29 ± 0.031 [17] (the analysis 
includes multipoles up to I = 400 in the W-band; a smaller evidence emerges from the V and Q bands). Secondly, it 
shifted the privileged direction very close to the ecliptic poles 0, [I?} ■ Although this near coincidence suggests that 
the asymmetry may not be cosmological, a systematical cause for it has not yet been determined. Specifically, ref. 
fl7| ruled out that the effect may be due to asymmetric beams, misestimated noise, or Zodiacal lights. 

In absence of a definite answer, it will be important to see whether the asymmetry will be also present in the 
forthcoming Planck data. In the meantime, it is also interesting to study whether there exist simple inflationary 
models which can (i) account for the asymmetry, and, possibly, (ii) provide other measurable predictions. As a first 



1 Ref. I13I studied instead the power spectrum for violation of translational invariance. A general framework for directional dependence 
of the power spectrum was studied in , where it was forecast that a quadrupole modulation of the power spectrum as small as 2% 
of the total anisotropic signal can be detected by the Planck satellite. 
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step towards this, ref. [18[ provided the formalism for computing cosmological perturbations on backgrounds with 
a privileged direction (Bianchi-I background, with a residual 2d symmetry) 2 . The formalism was then applied to 
the simplest case of inflation driven by a slow rolling scalar field, assuming that the expansion rate of one direction 
is different from that of the other two as an initial condition at the onset of inflation. One of the gravity waves 
polarization experiences a large growth during the initial anisotropic era (this is intimately related to the instability 
of Kasner spaces), which may result in a large B signal in the CMB [21]. Inflation however rapidly removes the 
background anisotropy. The modes that leave the horizon well after the universe has isotropized were deep inside 
the horizon while the universe was anisotropic, and one recovers a standard power spectrum at the corresponding 
scales. The signature of the earlier anisotropic stage are therefore present only on relatively large scales, which were 
comparable to the horizon when the universe was anisotropic. These scales can be visible today only if the duration 
of inflation is limited to a minimal amount. 

To avoid this tuning, one can obtain prolonged anisotropic inflationary solutions by introducing some ingredients 
that violate the premises of Wald's theorem [22j on the rapid isotropization of Bianchi universes. This has been realized 
through the addition of quadratic curvature invariants to the gravity action (23|, with the use of the Kalb-Ramond 
axion 24] 3 , or of vector fields 26]. A number of recent works focuses on this last possibility, which is probably 
technically simpler than the other two. 4 The underlying idea is that the vector has a nonvanishing expectation 
value along one spatial direction, causing that direction to expand differently from the other two. 5 In the standard 
case, i.e. for a minimal kinetic term (— F 2 /A) and no potential term for the vector, the vector vev rapidly decreases 
with the expansion of the universe, leading to a rapid isotropization. Therefore, the action for the vector needs to 
be carefully arranged. To our knowledge, four different possibilities have been explored. The first three possibilities 
are characterized by (i) a potential V (A 2 ) for the vector [26], (ii) a fixed spatial norm of the vector, enforced by a 
lagrange multiplier [10] , or (iii) a nonminimal coupling of the vector to the scalar curvature [30[ ■ 6 All these proposals 
break the U(l) symmetry which is present for a minimal vector action. This introduces an additional degree of 
freedom (the longitudinal vector mode) , which, for all of these models, turns out to be a ghost [ill, G2, H2 • This leads 
to instabilities of these models both at the linear (the linearized solutions for the perturbations diverge at a finite 
time close to horizon crossing) and nonlinear level (vacuum decay, with a UV divergent rate, into ghost-nonghost 
excitations). 7 

A completely different model was proposed in (36[. It is characterized by a scalar inflaton field, with a flat potential, 
and by a vector whose kinetic term is multiplied by a function of the scalar, — / (</>) 2 F 2 /4 . 8 As shown in 36], for a 
suitable choice of / the vev of the vector evolves slowly during inflation, and the model therefore supports a prolonged 
anisotropic stage. 9 Moreover, since the model is U(l) invariant, the problematic longitudinal vector is absent. 
One therefore expects this model to be stable. This was indeed shown to be the case in ref. [4l[. 10 We believe 
that this makes the model of [361 ] particularly interesting, since it is a complete counter example to the Wald's no 
hair theorem [22j, which is proven to be free of instabilities. In the present work, we study the phenomenology of 
this model. We compute the full spectrum of perturbations, and the resulting two point correlation functions. We 
particularly focus on two phenomenological signatures: the angular dependence of the scalar-scalar correlator, to see 
whether it can reproduce the value for reported in [l7j . and the scalar-tensor correlator, which is a distinctive 
prediction for this class of models (since it vanishes in the standard case). While, to our knowledge, the one obtained 
here is the first definite prediction for an anisotropic inflationary evolution, we show that the specific model proposed 
in [36[ does not reproduce the observed signal. 

The plan of this paper is the following. In Section [TT] we review the model of (36j and we study the background 



2 Ref. [19| provided the computation f or g eneral Bianchi-I backgrounds; the results of [18|| and [ill agree in the limit of 2d isotropy, and 
reduce to the standard computation [20| in the limit of 3d isotropy. 

3 For later studies of inflationary models with p— forms, see [25(1 . 

4 Several recent works [2^1 also discuss the possible non-gaussian signature from vector fields during inflation (with or without anisotropic 
expansion), as this can be a way to differentiate them from scalar fields. See also [28l | for an earlier study of the cosmological curvature 
perturbations generated by a vector field. 

5 While we focus on primordial inflation, vector fields with nonvanishing spatial vev have also been employed as sources of the late time 
time acceleration 1 2911 . 

6 Vectors with nonminimal coupling to the curvature on an isotropic inflationary background were studied in [3~D| . 

7 Refs. [33l l34fl ar gue d ag ainst this instability for models of the type (iii). We however believe that the arguments and the calculations 
presented in 111! |12|, |32|| are robust. For other stability studies, see [35|l . 

8 Ref. |37|| studied the generation of primordial magnetic fields in the case in which the vector is the electromagnetic potential. An 
analogous study was performed in |38H . 

9 A similar idea is also proposed in [3911 . where the function multiplying the kinetic term is taken to be an external function of time. For 
other works with vector fields with nonstandard kinetic terms during inflation, see [4d| . 

Ref. [4l| also presented the power spectrum of a subset of the perturbations of the model, including one of the gravity waves polarization, 
for which the computation is technically simpler than the complete one given here. 
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evolution. In Section |TT] we instead perform the computation of the perturbations and we show the resulting power 
spectra. The two main steps are the computation of the quadratic action for the perturbations (from which we obtain 
the linearized evolution equations for the modes), and the quantization of that action (which is needed for computing 
the initial conditions for the modes, and the correlators). For the latter, we employ the formalism of [44| for the 
quantization of a system of coupled bosonic fields. Our results are discussed in the concluding Section HVl 



II. THE MODEL AND THE BACKGROUND EVOLUTION 



We study the background solution in the model of [36[ , which is characterized by the action 

1 



/ 



d; y . 



-f) 



Ml 



V{4>) 



f(<t>) 2 F 2 



(2) 



Namely, there is a scalar field cj) with potential V ((/)), which is taken sufficiently flat to allow for a slow evolution of the 
scalar. There is also a vector field, which enters only with its kinetic term (F^ — d il A v — d v A^), so that the action is 
U(l) invariant. This kinetic term is multiplied by a function of the scalar, which causes a nonstandard evolution for 
the vector vev. For an appropriate choice of / (0) the energy density of the vector also evolves slowly, and the model 
admits a prolonged anisotropic inflationary solution (36j . 

We assume a Bianchi-I background with a residual isotropy in the plane perpendicular to the vector vev, 



ds 2 = -di 2 + a (t) 2 dx 2 + b {tf [dy 2 
and we parametrize the two scale factors by 



dz 2 



(AJ = (0, A 1 (t) , 0, 0) 



(3) 



(so that e Q = a 1 / 3 ^ 2 / 3 , e a = 6 1 / 3 a -1 / 3 ). Namely a parametrizes the overall volume expansion, while a controls 
the amount of anisotropy. More accurately, the normalization of the two scale factor is unphysical, and one can shift 
the values of a and a by a constant factor without changing the physics of the system. Therefore, the "initial" (i.e., 
at the initial time of a simulation) values ai n and a- m can be set to any arbitrary value. The degree of anisotropy is 
related to a /a: a constant a corresponds to a FRW flat geometry, with the Hubble rate given by d. 
The background equation of motion for the vector vev is 



A& + 2 



Ai=0 



(4) 



where dot denotes differentiation with respect to time, while /' - and, later, V' - denotes the derivative of that function 
with respect to the scalar field. This equation is solved by 



A x = p A e 



-a— 4cr 



(5) 



where pa is constant. We insert this solution in the other nontrivial background equations following from @. These 
equations then rewrite 



where we have defined 



3d 2 



3a 2 = 



1 

Ml 



v(4>) 



2a + 3d 2 + 3& 2 



3d a 



M 2 



3M p 2 /(0r 



4> + 3d j> + V ((f>) = 



/w 3 



2 /(</>)' 

v (<!>)- 



6/(</>r 



(6) 



PA{t) 



e- 2a ^ PA 



(7) 



4 



The first of ([6]) is the tt Einstein equation for the system, and, for a = Pa — 0, it reduces to the standard Friedmann 
equation. The last term in this equation is the energy density of the vector field. Since the normalization of the scale 
factors is unphysical, the quantity Pa must be independent of it. Therefore, if one wishes to change the "initial" 
values of the two scale factors, the integration constant pa must also be changed according to pa oc e 2ain+2<7in . rj^g 
second and fourth equation in (|6]) are a combination of the spatial xx and yy = zz Einstein equations, while the 
fourth equation is the equation for the scalar field. One of these three equations can be obtained from the other two, 
and from the tt Einstein equation, as a consequence of a nontrivial Bianchi identity. 

We are interested in an inflationary background solution of © in a regime of slow roll and small anisotropy. The 
first and last of © can be approximated as in the standard case, and combined to give a (4>) ~ — M Y V , dip. We 

can also find how / and V need to be related to each other to have a prolonged anisotropic stage. Namely, we require 
that the ratio between the energy densities of the vector and the scalar field remains approximately constant during 
inflation. In the slow roll regime, and for small anisotropy, we have 



9 — An 

PA _ P A C 

P4> ~ 2V(4>) V/WO 



(8) 



and, given that the first factor is very slowly evolving in this regime, we require [36j that / w exp [—2a] 



exp 



For definiteness, we will consider the simplest chaotic inflationary potential 

F(0) = im 2 2 , /(0)= eX p(^O (9) 



2" T 7 J VT ' \2M 2 y 

where c is a numerical constant close to one. 

We expect that, for c > 1, / decreases more quickly than e~ 2a during inflation, so that the anisotropy actually 
increases; to verify this, we study the first and last of ([6]) in more details. Namely, we disregard the terms proportional 
to a and Pa in the first equation, and the term proportional to <f> in the last one (it is consistent to keep the term 
proportional to Pa in the last equation, and disregard it in the first one, provided that <f> MjJy/c, which is indeed 
a good approximation during inflation). These two equations can be then integrated to give [361 ] 

e c 4 +4Q « J- -||^ + D e -4(c-i)a (10) 

where D is an integration constant. For c < 1, the second term on the right hand side of this equation increases over 
the first one during inflation. If this term is dominant, then Eq. (fTO)) reduces to a s» — <fr 2 / (4 M 2 ) ( plus an unphysical 
constant), which is the standard isotropic result. On the contrary, for c > 1, the first term of (Il0|) increases over the 
second term over time. If this term is dominant, Eq. (p~0|) reduces to a w — c0 2 / (4Af 2 ). Combining this with the 
approximate tt Einstein equation, 6d 2 ~ m 2 2 /M 2 , gives 1365 

m 2 6 , , , 

3a0w (c> 1 (11 

c 

which gives a value for (j) about 1/c times the standard result. To quantify the anisotropy, we neglect the a term in 
the third of Combining the resulting expression with (|10p . in a regime in which the term proportional to D can 
be disregarded, we find 



a 
a 



2 c - 1 M 2 



(c > 1) (12) 



which indeed confirms that the anisotropy increases during inflation. 

This quantity is approximately equal the ratio between the energy densities of the vector and the scalar. Indeed, 
combining eqs. © and (ITU1) . in a regime in which the term proportional to D can be disregarded, we also find 

PA ^ C-l M 2 . \ / 1 r,\ 

— ^—^--nr ( c>1 ) ( 13 ) 

Namely, for c < 1 the system evolves towards isotropy, while for c > 1 the system evolves towards isotropy during 
inflation. In both cases the solution is an attractive one. One can decide to take c < 1 and start away from the 
isotropic attractor solution. If c is sufficiently close to one, the anisotropy will decrease very slowly, and still give some 
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observable nonstandard signature. The result will however be sensitive on the initial conditions, and not only on the 
model. On the other hand, for c > 1 we can start in the attractor anisotropic solution characterized by (fTTj). The 
underlying idea is that inflation lasted much more than the observable last 60 e-folds, and that the solution converged 
to the attractor one during that time. In this case, the phenomenological signatures of the model are insensitive on 
the initial conditions, precisely as in the standard inflationary case. 11 For this reason, we only study the c > I case 
in this work. 

We conclude this Section with two remarks. Firstly, we note that the anisotropy is proportional to c — 1. The 
anisotropic attractor solution is continuously connected to a FRW solution in the c — > I limit. We expect standard 
results for the perturbations in this limit, as the computations presented in the next Sections confirm. Secondly, while 
the anisotropy increases during inflation, it decreases after inflation. Indeed, after inflation <f> oscillates around zero, 
with a decreasing amplitude. Then / —¥ I, and the mechanism of prolonged anisotropy becomes ineffective. The 
amplitude of the vector rapidly decreases, and the background evolution becomes isotropic. 

For illustrative purposes, we show in Figure Q] the evolution of the anisotropy a /a as a function of the number of 
c-folds N = a, normalized to zero at the end of inflation. 12 We show the evolution for two different values of c, 
starting from the slow roll anisotropic initial condition. 13 We also show the evolution of the anisotropy parameter 
given by the slow roll solution, eq. (JUJ). We see that the slow roll expression is very accurate during most of the 
inflationary evolution. 



III. PERTURBATIONS 

This Section studies the perturbations of the model [36[ discussed in the previous Section. The discussion is 
divided in several Subsections. In Subsection IIII Al we introduce the perturbations, and we classify them according 
to how they behave with respect to 2d spatial rotations in the isotropic yz plane. We also perform the gauge choices 



Even in the standard case, one can assume that the inflaton was not yet in the attractor solution when the largest observed multipoles 
left the horizon. For instance, a fast roll evolution at that stage results in a suppression of the CMB quadrupole [42| . This signal is 
however dependent on the assumed initial conditions. 

Since a <g a, we define the number of e-folds N and the end of inflation as in the FRW case. This gives N = a, and the end of inflation 
occurs when a + a 2 becomes negative. 

Specifically, we use the last three equations of l[6j in our numerical evolutions. We also satisfy the first of $6% through the initial 
conditions (if this equation holds at the initial time, it is preserved by the remaining equations): we relate the initial value of <j> and of 
<t to that of a through eqs. mil l and (112b ; we insert these expressions in the first of ((BJ, and we then use this equation to relate di n 
to the initial value of 4>. In this way, all the initial conditions are given in terms of </>; n , which completely specifies any point along the 
attractor solution. 
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which completely fix the freedom associated with general coordinate transformations, and with the U(l) invariance 
of the vector. We then integrate out the nondynamical perturbations, and we provide the quadratic action for the 
dynamical modes. In Subsection IIIIBI we quantize this action, and we define the initial adiabatic vacuum, valid 
during inflation in the sub-horizon regime. In Subsection 1111 CI we introduce the two point correlation functions for 
anisotropic backgrounds, and we express them in terms of angular dependent power spectra. In Subsection 1111 Dl we 
write more explicitly the initial conditions for the modes (coming from the adiabatic vacuum obtained in IIIIBj) and 
their evolution equations (coming from the action obtained in IIII A[) . Finally, in Subsection IIII El we discuss which 
combination of perturbations reduce to the standard one as the universe isotropizes, and we provide their power 
spectra. 

Both the computation of the quadratic action, and its quantization, follow formalisms that have been developed 
elsewhere. Although we made an effort to keep the present discussion self contained, some details have been omitted 
here for brevity. The interested reader can refer to appendix B of and to Section III of [HJ for more details on 
the computation of the quadratic action, and to Section II of (iij for the quantization of the perturbations based on 
this action. 



A. Classification, gauge choice, and quadratic action for the dynamical modes 

The most general set of perturbations about the anisotropic background solution discussed in the previous Section 
is given by 

/-2$ ad xX bidiB + Bi) \ 

%„= -2a 2 * abd x (diB + B^ , SA„ = (SA , SA lt drfA + SAi) , H (14) 

V b 2 [-2Z» tf + 2E& + E (itj) ] J 

where the index i — 1,2 spans the isotropic yz plane, and -E(ij) = diEj + djEi. With this choice, the perturbations 
are decomposed in scalar and vector with respect to rotations in the isotropic plane. This is a convenient procedure 
since the two different types of modes are decoupled from each other at the linearized level, and can be studied 
separately [HI]. 14 Specifically, the modes Bi, BiE%, SAi are 2d vector modes, satisfying diBi = • • • = 0, and they 
encode one degree of freedom each. The remaining modes are 2d scalars, which also encode one degree of freedoms 
each. Altogether, we have 15 degrees of freedom in the perturbations (fT4|) . 
We Fourier transform the perturbations, 

/rl 3 k 
-—e^5(t,k) (15) 
{2ny /2 

where 8 denotes any of the perturbations, and we use the same symbol for a perturbation in real and in momentum 
space. The reality of 8 (x) imposes the condition 8' (k) = 8 (— k) (we use hermitian conjugate, and not simply charge 
conjugate, since in the next Section the modes are treated as operators, in order to quantize them). At the linearized 
level, all modes are decoupled from each other, and can be studied separately. We denote the comoving momentum of 
the mode as k = (fcj,, kx2, fcra)- Due to the symmetry in the y — z plane, physical results depend on kx = \J k T2 + k T3 
rather than on kx2 and kx3 separately. 15 We denote the physical momentum of the mode by p = {pl, PT2, Pts), 
where pl — k^j a and pxi = kxi/b. As for the comoving momentum, we define px = \/Pt2 Pt3- 

It is convenient to write write explicitly the single degree of freedom encoded in the 2d vector modes. In momentum 
space, we have 

II, i ',,/. ; , B v , Bi = ieijkrjB v , E l = id ld k T jE v , SAi = ieijk T j SA V (16) 
where is antisymmetric, and ei2 = 1. 



This strategy is mutilated from the one adopted in the standard 3 dimensional case, in which the modes are classified in scalar, vector 
and tensors with respect to 3d rotations. Two dimensional tensor modes do not exist, since imposing the transversality and traceless 
conditions eliminates all the degrees of freedom of a 2 X 2 symmetric tensor. 

This was exploited in [Till , where kxs was set to zero. As a consequence, the perturbations in real space do not depend on z, and the 
decomposition i 1 41) was simpler. Here we adopt the more general formulation used in |2l|. where both kx,2 and fcy,3 can De nonvanishing. 
We also point out that spatial derivatives are introduced in the parametrization (114 I t for algebraic convenience. In practice, we assume 
that both and fcy are nonvanishing. This restriction does not affect the computation of any observable, since modes with = 
and/or kj- = constitute a subset of zero measure when one integrates over the momentum to obtain results in real space. 
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To proceed, we need to fix the gauge freedoms of the system. We start from the freedom associated with general 
coordinate transformations. To do so, we can either choose a gauge that completely removes this freedom, as done in 
[3 or we can rewrite the action and the equations for the perturbations in terms of gauge invariant modes, as 
done in [12, H^]. The two procedures are equivalent. In the present work, we choose the first one, which is algebraically 
simpler. Specifically, we set 



7li,2ds 







(17) 



which, in the parametrization (fT4)) . gives B = £ = E = Ei = 0. One can check (see Appendix B.l of [18|) that 
indeed (i) this choices can be always made, and (ii) it fixes completely the freedom associated with thee coordinate 
transformations. There is also the U(l) gauge associated with the transformations Ap — > + dp£. We fix this by 
setting 8 A = 0. 

This gauge fixing leaves us with seven 2ds modes ($, x, B, ^, 5Ao, SAi, 8<f)) and three 2dv modes (Bi, Bi, SAi). 
Not all these modes correspond to physically propagating degrees of freedom. Indeed the modes 5A and Sgo^ enter 
in the quadratic action of the perturbations without time derivatives [ill]. As a consequence, the equations of motion 
for the perturbations are algebraic in them, and their value is specified in terms of the values of the other modes, 
without introducing independent degrees of freedom. Our gauge choice is motivated by the fact that it preserves the 
SAo and dgoii perturbations, so that the identification of the nondynamical modes is immediate (in other gauges, the 
nondynamical modes correspond to more complicated linear combinations of the perturbations which are preserved 
in those gauges). The nondynamical modes need to be integrated out of the action. 16 Namely, we express them in 
terms of the dynamical modes (through the corresponding Einstein equations), and we insert these expressions back 
into the action. In this way, we are left with an action in terms of the dynamical modes only 17 



7 (2) 

pcrts. 



(2) 



2dv 



The fields entering in these actions are not canonically normalized, 
through the redefinitions 



B h 5A t 



(18) 



The canonically normalized fields are obtained 



Oil 



e 2 



in the 2ds sector, and 



V2M p (a + a 



a 2p 2 d + (2p\ - y\) & 
V2M p p^ (a + a) 



PA 



f{4>)PT y/2M p f(<j>) 2 (a + a) 



-H, 



(19) 



Bi 



-V2eijkTje ^ 



-13 J 



M p k L k\ 



in the 2dv sector (where ey is antisymmetric, and eia = 1). 

Once expressed in terms of these fields, the two actions in (fT8|) rewrite 



(20) 



3(2) 
7 2ds 

3(2) 
'Idv 



| / dtd 3 k 



dtdfk 



)\ Y. + Y l\A: - Y}K S Y S - ^ ()2 Ys 
Y^Y + y ^ K Y —Y^KY — Y^ Q 2 Y 



(21) 



This is what is also done in the standard computations [20|| . For instance, in the case of a single scalar inflaton on a FRW background, 
the number of perturbations is 11 (initial perturbations in the metric and in the scalar) —4 (after gauge fixing) —4 (after eliminating 
the nondynamical modes) = 3 (namely, the scalar density contrast, and the two gravity waves polarization). Our procedure extends this 
computation to the more general background we are studying. In Appendix [B] we verify that our computation reduces to the standard 
one in the limit of isotropic background. 

The procedure of integrating out the nondynamical modes is described in details in Section III of [32H . 
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where 




y ( f> ; (-2. 



and where the explicit expressions for the matrices K S:V and ft 2 v are given in Appendix lAl 



B. Quantization and initial adiabatic vacuum 



We need to quantize the two actions (|2Tj) in order to provide the initial conditions for the modes and the expressions 
for the correlators. We do not need to discuss the two systems separately, since the two actions are formally the same, 
and the fields and the matrices entering in them have identical properties. In both cases, the actions are of the type 



dtd A k 



y^y + yt ky - y+ KY -Y^Q 2 Y 



(23) 



where Y is an array of fields, K a real and anti-symmetric matrix, and fi 2 a real and symmetric matrix. These 
matrices are unchanged under the parity transformation k — > — k (as can be seen from the explicit expressions 
given in Appendix |A"]) . It can be checked from the reality condition stated after eq. (|15p . and from the definitions 
(fT9| and (|20| of the canonical modes, that, for both systems, any of the Yi fields entering in the array Y satisfies 

y/(k) = y(-k). _ 

To remove the mixed terms proportional to the matrix K , we first perform the field redefinition 



ip = RY 

where R is an orthogonal matrix (so that Y^Y = ip^ip ), satisfying 

R=RK , i?latc = l 



(24) 



(25) 



where the second condition states that R should reduce to the identity at late times, when the universe becomes 
isotropic, and K — > (in which case, the rotation (|2"4"|) is no longer needed). The first condition can be also written 
as K = R T R; we also note that R (— k) = R (k), since this property is also satisfied by K. As a consequence, each of 
the fields entering in the array ip satisfies ip\ (k) = ipi (— k) . In terms of the fields ip, the action (|23|) rewrites 



S 



dtd 3 k 



tr>2, 



n 2 = R(n 2 + K T K) R T 



(26) 



where we note that Cl 2 is real, symmetric, and invariant under k — > — k. 

We can also define the real space fields if> (t, x) as in ([T5|) . These fields are real, and their action is formally identical 
to the action for the coupled bosonic system quantized in 44]. Therefore, we quantize the fields ipi as done in that 
work. 

We introduce the matrix C satisfying 



C 1 W C = diag [u(, 

Ccnd = 1 



) 



UJ 



(27) 



where the second condition follows from the fact that, in the late time isotropic limit ft 2 = il 2 is already diagonal. 
We note that C is orthogonal, and unchanged under k —> — k. We then define [3] 



fl>i (k) 
m (k) = ip\ (k) 



h jt (k) m (k) + h* jt (k) a] (-k) 
hfl (k) a, (k) + hi (k) a] (-k) 



where a and a) are annihilation and creation operators, respectively, satisfying 

(k) , a] (k')l = S 3 (k - k') Sij 



(28) 



(29) 
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From the equations of motion following from (l26l) . and from the fact that iri = ipi, we find that the coefficients hij 
and , obey the evolution equations (in matrix from) 

k = h-Th , h= -rh-cu 2 h , T = C T C (30) 

From the parity properties of the matrices C and lo, and from the initial conditions (which we determine below, see 
eq. (|3"5)) ). we see that and are unchanged under k — > — k. 18 
We further define 

h=-±= (a + p) , h==£L{tx-P) (31) 



It has been shown in [44j | that the normal ordered hamiltonian for the fields ipi can be then cast in the form 

H = f SkLo.bl (k) 6(k) (32) 



where hi and b\ are new annihilation and creation operators, related to those defined in (|28p by (notice that the 
matrices a and /3 are unchanged under k — > — k) 



S(t,k) \ = / a /3* \ / a(k) 
S+(t,-k);-^ ^ ) ^fit(-k) 



(33) 



We see that the hamiltonian is diagonal in the 6i,&] basis, so that these operators annihilate and create quanta of 
the (time-dependent) physical eigenstates of the system. The matrices a and (3 generalize to a system of N coupled 
fields the Bogolyubov coefficients that are needed for the quantization of a field with time dependent frequency. As 
shown in [44j |. the canonical quantization of the ipi fields imposes the conditions 



aa 



f3*/3 T = 1 , a/3 t - f3*a T = (34) 

Moreover, from the evolution equations (IBTfl) . and the definitions (|3"Tj) . one finds that a and j3 obey the evolution 
equations [44| 

a = —iuioi + 7^-/3 — la — J/3 

(35) 

/? = ioj(3 + £-a - 1(3 - Ja 



/ = ± ( V^r_L + _Liy^) J = l{ y /zr-L--Lr y fr) (36) 



where 



An inspection of the initial matrices K and fl 2 shows that, at early times (when the mode is deeply inside the 
horizon) f2 2 = p 2 1 + O (H) and K — O C\/c — 1 H) (in these expressions, H = a/a ~ 6/6 under the assumption of 
small anisotropy that we are making in this work). As a consequence, 

uj 2 ~p 2 + {H 2 ) , I\ I, J, - = O (J?) (37) 

w 

in this early time regime. Therefore, we can disregard all but the first term in both of the right hand sides of (|35[) . 
This leads us to the initial adiabatic vacuum solutions 19 

Nearly = e" 4 dt " , /3 oarly = (38) 

This solution obeys (|35[) in the early time regime, satisfies the quantization conditions (|34[) . and corresponds to an 
initial empty vacuum, if we impose that the vacuum state is annihilated by the operators entering in 



This is why we wrote h*- (k) and h*- (k), rather than h*. (— k) and h^ - (— k), in the decompositions 1281 1. 

Notice that « car i y is diagonal; the allowed initial conditions are actually more general than 1381 1. since one can multiply each diagonal 
entry of a e arly by a constant, and arbitrary, phase factor e I7i . This amounts in changing the matrices a and given here by a matrix 
multiplication from the right, a — ¥ aP, — > P, where P = diag (e' 71 , . . .e l7w ). The equations of motion 1351 1 are unchanged by 
this multiplication. The same is true in terms of the matrices h and h. Since h enters in the observable two point correlation function 
through the combination hh', see eq. 1421 1. the matrix P drops from the observable result. This confirms that the arbitrary phases 
contained in P are unphysical, and can be set to any value. We use this freedom to set all the phases to zero at the initial time of our 
numerical simulations. 
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C. Two point correlation functions 

By combining the various redefinitions given in the above Section, we can write the canonically normalized fields 
in real space as 

— p e* kx [T« (t, k) a, (k) + T*j (t, k) a] (-k)J (39) 
T^(k)= (R T Ch).. (40) 



where 



(notice that h and T coincide at late times). 

The (statistically averaged) two point correlation function can be expressed as the quantum expectation value 

dj (x, y) = ±(Yi (t, x) Yj (t, y) + Yj (t, y) F 4 (t, x)) (41) 

where the symmetrization is required since the statistical average is a classical operation, independent of the ordering 
chosen (when computing any correlation with real data, Cij (x, y) = Cji (y, x)). 

We insert (|3U)) into (|4ip. The resulting expression can be then simplified using the commutation relations (j2"9")l and 
the fact that the vacuum state is annihilated by hi at all times (since the quantization of the previous Subsection is 
performed in the Heisenberg picture). After some algebra, we find 

M*,y) = / ^3e^-)Re[(TTt)J (42) 

We now define the power spectra associated with these correlators. All of them are of the type 

c -=/!r' Ml -^ (k) (43) 

where the function T is real and it depends only on the absolute values of the components of k along the anisotropic 
a;— direction (which we denoted by |A;j,|), and on the y — z plane (denoted by kx — y/kj, 2 + k^ 3 ). Thanks to this 
property (which follows from the symmetry of the background under rotations in the yz plane, and under parity) , for 
any two points x and y, we can always choose the y and z axes of the system such that the third component of x — y 
vanishes (without changing F). We therefore set 

r = x y = (r L , r T , 0) (44) 

and k = k \f\ — £ 2 cos (f>k, y/l — £ 2 smfikj in the integral (|4"3"]l . where £ is the cosine of the angle between the 
x— axis and k. The function J- does not depend on and is even in £. The integral over then gives 

C? = J® J d£coa(ktr L )J (ky/T=?rT) Pr (45) 
where J is the Bessel function of the first kind, and where we have introduced the power spectrum 



^=^2^(*.0 ( 46 ) 

27T Z 



In the case at hand, the power spectra depend both on the magnitude of the momentum of the modes, and on the 
angle between the momentum and the anisotropic direction. 

On anisotropic backgrounds, the power spectrum is isotropic, and eq. (1431) reduces to the standard expression 

f dk sin (kr) ^ , 
Cjr= — —^r 1 V ? isotropy (47) 
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D. Evolution of the perturbations and initial conditions 



The equations of motion for the dynamical perturbations follow from (|2ip . We expressed the action in momentum 
space, by Fourier transforming the starting modes as in (|15p . and by introducing the canonically normalized fields in 
(TH)1) and (j2"0")l . The equations of motion for the coefficient of the Fourier transforms of the canonical fields are 



Y.+2K.Y.+ lSr.+K.)Y.=0 



nl + k v )y v = o 



(48) 



where the matrices K s v and tt 2 v are given in Appendix [A"l 

We stress that these equations are a closed subset of the linearized Einstein equations for the perturbations. Specifi- 
cally, all the perturbations of the model can be divided in dynamical and nondynamical ones. The nondynamical ones 
enter in the Einstein equations without time derivatives. One can solve the Einstein equations for these perturbations, 
and express the latter in terms of the dynamical ones. One then insert these expressions into the remaining Einstein 
equations. The resulting expressions coincide with eqs. (|48[) . 20 

We need to exactly specify what the coefficients Y s ^ and Y v ^ exactly are. The standard way to compute the 
generation of perturbations dur ing inflation is a semiclassical computation, in which the perturbations are quantum 
fields on a classical background [20]. We take this approach in this work. In most systems, the canonical perturbations 
are decoupled from each other, and can be quantized separately. This is not the case for the system we are studying, 
and the quantization had to be done accordingly. In Subsection IIII Bl we introduced an array of annihilation/creation 
operators, and we saw that the Fourier coefficients of (|15l) are actually linear combinations of these operators, cf. eq. 



(49) 



Y (t, k) = Tij (t, k) aj (k) + (t, k) a] (-k) 
Inserting this decomposition into (|48p. we find 



2K il T lj +lW + K) T 



-ij 



Y*j + 2KuT*j + IQ + K) T 







(50) 



both in the 2d scalar and 2d vector sector. The linear combinations multiplying different annihilation and creation 
operators need to cancel separately; therefore 



+ 2K a Tij 



T 



ij 







(51) 



These equations also guarantee that the linear combinations multiplying the annihilation operators vanish, since the 
matrices K and f2 2 are real. 

Having determined the evolution equations obeyed by , we now turn to the determination of their initial condi- 
tion. We find 



T = R T Ch , f = -KR T Ch + R T Ch 



(52) 



The first expression is simply the definition (|40[) . The second expression is obtained by differentiating the first one, 
and by using the first and third of (|5U|) . as well as K = —R T R . 

The initial values of R T C and w are obtained by the diagonalization of ft 2 + K T K at the initial time. Indeed R T C 
is defined as the matrix that diagonalizes il 2 + K T K, while to 2 is the diagonal matrix formed by the eigenvalues of 
n 2 + K T K, see eqs. (|2"o| and (|2"T|) . The initial values of h and h follow instead from eqs. (f5T|) . and from the initial 
conditions for a and f3 according to the adiabatic vacuum prescription, cq. (1381) : a- m = 1, (5- ln = 0. 



E. Power spectra after isotropization 

As shown in Appendix |B] as the universe becomes isotropic after inflation, the canonical perturbations that we 
have introduced in IIII Al become the standard scalar and tensor modes of FRW cosmology. More precisely, we find 
that, in this regime, 

r = ^j: v + ' h + = -^rr H + > h * = -&r!or H x ( 53 ) 

a? 1 2 4> a A/ M p a*' 1 M p 



The explicit proof of this is given in Section III of 13211 for details. 
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where 1Z is the standard (scalar) comoving curvature perturbation, and h + and h x are the two standard gravity wave 
polarizations. 

We also shown in Appendix [B] that, in the isotropic regime, our formalism reproduces the standard evolution 
equations for these modes. These equations are decoupled, so that the modes evolve independently from each other 
(at the linearized level). However, the two 2d scalar modes V+ and H + are coupled to each other, and to the 
mode A+ during inflation, when the background is anisotropic. The mode H x is always decoupled from these three 
modes, but it is coupled to the mode A x during inflation. The coupling modifies the diagonal correlation functions 
(V+), (A+), (H x ), (A^) with respect to the standard inflationary results, and it introduces the nondiagonal 

correlations (V+ H+), (V+ A + ), (H + A + ), (H X A X ), which are absent in standard inflation. 

As the universe isotropizes, the two modes A + and A x become the two transverse polarizations of the vector field. 
These modes rapidly decrease after inflation (when the evolution of the vector becomes standard). Therefore, all 
correlators involving these modes become negligible at late times, and we disregard them in the remainder of this 
work. The remaining correlators can be written as in eq. (|45[) . in terms of the power spectra 



P ™ = if f 3 (T,Tt) u 

= ^^M^ (54) 

To evaluate the power spectra, we impose the initial conditions on T as discussed at the end of the previous Sub- 
section. We then evolve T through their equations of motion (fSTj) . We remark that the resulting power spectra are 
dimensionless, and are written in terms of only physical quantities (namely, they are insensitive to the normalization 
of the scale factors; this is because only physical quantities appear in the initial conditions and the evolution equations 
for T). 

Our results are shown in Figures [2] [3l and [4] Figure [2] shows the evolution of the power for the specific value of 
c — 1 = 10~ 5 in the kinetic function ©, and for a specific mode: we denote by k* the magnitude of the momentum 
of a mode which barely exits the horizon at the end of inflation; we choose a mode with momentum k = 10~ 20 fc*, 
and we choose the cosine of the angle between the momentum and the privileged direction to be £ = 1/2 (such values 
have no particular meaning, and are just chosen for illustrative purposes; different values of k and £ lead to the same 
qualitative behavior); this mode is initially deeply inside the horizon, and it leaves the horizon about 50 e- folds before 
the end of inflation. 21 The power in the scalar-scalar and tensor-tensor correlations behaves analogously to the 
isotropic case. The power in the scalar-tensor cross correlation is instead very small initially, and slightly increases 
during and immediately after inflation. The power becomes constant (at the value seen in the latest time shown in 
the Figure) as the universe isotropizes after inflation. 

In Figure|3]we show the power spectra for the specific value of c— 1 = 10~ 5 in the kinetic function ©. The smallest 
momentum shown corresponds to modes that exited the horizon about 60 e-folds before the end of inflation. For 
the scalar-scalar and tensor-tensor case, the standard result is also shown for comparison. The scalar-scalar power 
spectrum (bottom left panel) is slightly greater than in the isotropic case; this can be compensated by decreasing 
the scalar field mass (for this reason, the ratio m/M p has been kept as a free parameter in the Figure). The angular 
dependence of Pun is of O (10 _1 ) at the largest scales, while it slowly decreases at greater scales. The tensor-tensor 
power spectra (bottom panels) are closer to the standard result, and they exhibit a much milder angular dependence 
(we found g* for the scalar spectra is suppressed with respect to in the tensor spectrum by approximately the 
ratio between the power of the tensor and the scalar spectra). Moreover, the results for the two polarizations are 
nearly identical. The scalar-tensor cross correlation (upper right panel) shows a stronger angular dependence, but it 
is smaller than the other two spectra. 

In Figure H] we show the angular dependence g* (defined in equation ([!}) for the scalar-scalar power spectrum, for 
different values of c in eq. ([SJ . We remark that all values of c provide a negative g* (in Figure H] we actually show 
|<7*| = — in logarithmic scale). The value of g* shown in the plots is obtained by comparing, for each value of 



21 Since this choice of c ~ 1 corresponds to a very small anisotropy, we define the horizon, and the number of e-folds, only through the 
evolution of the "average" scale factor e a , as in FRW cosmology. 
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10 -3 I , y , , , , , 1 

-50 -40 -30 -20 -10 

e-folds 

FIG. 2: Time evolution for the power of a specific mode, on a nearly isotropic background (c — 1 = 10~ 5 ). The number of 
e-folds a is used as a "time" variable, and it is normalized to at the end of inflation. The mode shown leaves the horizon 
about 50 e-folds before the end of inflation. See the main text for details. 



k = |k|, the power at £ = 0.1 and at £ = 0.9. We have however verified that the ACW parametrization ([T} is very 
accurate, in the sense that, once g* and P (k) are obtained from the results at £ = 0.1 and £ = 0.9, also other values 
of £ are fitted very well by ([1]) . 

The significance of the results summarized in these Figures is discussed in the next concluding Section. 



IV. DISCUSSION 

In this work we compute the precise phenomenological signatures for the model of [36| . This model admits an 
anisotropic inflationary background evolution, supported by the combined presence of a scalar and a vector field. 
This solution is (mathematically) continuously connected to an isotropic solution, in the limit in which the energy 
density associated to the vector is sent to zero. Moreover, the model is free of the ghost instabilities that plague 
other models with vector fields during inflation, due to the U(l) invariance of its action under a shift of the vector. 
Therefore, it offers a complete, and stable counterexample to Wald's no hair theorem on the isotropization of Bianchi 



spaces [22 1 . 



We studied the simplest realization of the idea of [36j, in which a single vector is present, and the potential of the 
scalar is taken to be that of massive chaotic inflation. We found that the scalar-scalar correlation function exhibits 
an angular dependence which is however of the wrong sign to account for the breaking of rotational invariance seen 



in the data: the model gives a negative value for g*, while the analysis of [17[ indicates that = 0.29 ± 0.031 
in the WMAP W-band. We also found that the amount of anisotropy in the spectrum (the order of magnitude of 
\g*\) is not of the same order of magnitude as the amount of anisotropy in the expansion (the order of magnitude 
of AH/H ^ a/a, where H is the average expansion rate, and AH the difference between the expansion rate in the 
different coordinates). A = O (l0 _1 ) is obtained when the anisotropy in the expansion is of O (10- 7 - 10~ 6 ) 
during inflation. It is also worth noting that the anisotropy in the tensor-tensor spectrum is much milder than that in 
the scalar-scalar spectrum (we find that the suppression is approximately proportional to the ratio between the two 
power spectra). 

A distinctive feature of the model (and, of anisotropic spaces in general) is a nonvanishing scalar-tensor correlation, 
which, if sufficiently high, may be detected through temperature-B mode correlation in the CMB. A naive estimate 
actually suggests that the scalar-tensor correlation could be higher than the tensor-tensor one. Indeed the amplitude 
of a tensor mode is approximately multiplied by 2 ^fi with respect to that of a scalar mode, where e = —H/H 2 is a 
slow roll parameter. 22 The tensor-tensor correlator is suppressed by one more power of this factor than the tensor- 



Since we are considering a small anisotropy, we can use the FRW computation in this estimate; then, comparing the two equations in 
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k/k. k/k. 



FIG. 3: Scalar-scalar, scalar-tensor, and tensor-tensor power spectra for a nearly isotropic background (c — 1 = 10~ 5 ). Notice 
that the scale on the y axis is different for the different panels. For the scalar-scalar, and the tensor-tensor spectra, we also 
show the standard result for comparison (in the present model, the isotropic limit is reached for c = 1). 

scalar correlator. However, the scalar-tensor correlator vanishes for an isotropic space, and hence must be suppressed 
by some factor related to the asymmetry. If one naively assumes that this factor is \g*\, one would get the prediction 
that Pnh/Phh — \g*\/ (2\/e)- O ne would then get Pnh/Phh — 6\g*\ for chaotic inflation, and a greater value for other 
inflationary models, characterized by a smaller value of e. Our explicit results show that this estimate is reasonable, 
but not exact. For the case of c— 1 = 1CP 5 shown in Figure [3] (giving ~ —0.23), the estimate gives Pnh/Phh — 1-4, 
while the actual spectra give Pnh/Phh — 0.05 — 0.25, depending on the orientation of the mode (the value of £). 

To conclude, while the simplest realization of [36| cannot explain the breaking of rotational invariance seen in the 
data, and, most likely, it does not give rise to an interesting scalar-tensor correlation, this model is, to our knowledge, 
the first complete and stable model of anisotropic inflation for which the phenomenological predictions strictly follow 
from the action (and not from arbitrary initial conditions), and have been computed. As the observed breaking of 
rotational invariance awaits for a confirmation, or a refutation, from Plank, our work provides the tools for studying 
different models, to see whether they can reproduce the WMAP feature, and perhaps lead to new predictions. 



JB9) we see that |V+| ~ \H +lX \. This gives \h +lX /H\ ~ y/2cj>/ (M p H) ~ 2e. 
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FIG. 4: Angular dependence (expressed through the ACW g„ factor) of the scalar-scalar power spectrum for different choices 
of c. 



Note added: After the completion of the analysis reported here, and as we were finalizing the preparation of 
the present manuscript for submission, ref. [45| appeared on the preprint archive, which also computed some of the 
correlation functions computed here. Why we obtained the expressions for the correlators analytically, and we then 
evaluated them numerically, ref. H performs a fully analytical computation, treating the anisotropy as a small 
perturbation. The relation between and a/a found in [45[ and in the present work are in very good agreement 
with each other. The two studies also agree on the fact that the tensor-tensor correlator shows a smaller angular 
dependence than the scalar-scalar one. The scalar-tensor correlator has not been computed in 45 1. 
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Appendix A: Explicit quadratic action of the perturbations 



We provide here the explicit expressions for the matrices entering in the action and the evolution equations of the 
canonical modes, eqs. (|2H and 
For the 2d scalar modes we have 



K„ = 









K 13 , s 
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K 13 , s 


-K 23 , s 






o 2 
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and where we have defined 

For the 2d vector modes we have instead 

K„ -- 



V=2p 2 a+ (2p| - p 2 T ) a 
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The physical momenta entering in these expressions are related to the comoving ones given in the main text by 
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Appendix B: Gauge invariant perturbations in terms of perturbations in our gauge, and late time isotropic 

limit 



In this Appendix we discuss the late time interpretation of the perturbations in the gauge chosen in the main text. 
This interpretation is done when the background has become isotropic, and the vev of the vector has gone to zero. 
In this case, the Fourier coefficients of our metric perturbations read 
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as can be seen by imposing the gauge (JTTJ) on the parametrization given by (fTI)) and (fTBjl. We recall that (kL, &T2, krs) 
denotes the comoving momentum of the mode we are studying. 23 We perform a rotation to a coordinate system for 

which the momentum is along the third direction, fc M — >• k^ = R 1 ^ k u — ^0, 0, \fi? L 

of the rotation matrix is 



~~k~T2~ 



The explicit form 
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(B2) 



In the main text, the suffices L and T refer to "longitudinal" or "transverse" with respect to direction x, which was the anisotropic one 
during inflation. 



18 



where kx = \J k^ 2 + kj, 3 , and the metric transforms as 

Qiiv -> SVtr/ = -R" 3a/3 -Rf 



(B3) 



Although the gauge in which eq. (jBip appears is nonstandard, we can combine the metric perturbations (|Blj) into 
the gauge invariant expressions that are commonly used. These gauge invariant combinations are usually expressed 
starting from the most general metric perturbations, classified as scalar, vector, or tensor, with respect to 3d spatial 
rotations: 



5g 00 = -2$* , 5g i = 2a(B„ i + d i B t ) , Sg i3 ; = a 2 [-2** S zj + 2 E* 



E, : 



(B4) 



where i = 1, 2, 3. The (3d) vectors B„i and E*i are transverse, while the tensor mode h*ij is transverse and traceless. 
The remaining modes are scalar, and are coupled also to the perturbation of the scalar field <50*. Out of these modes, 
we are interested in the the gauge invariant scalar combination [20L |43j 



(B5) 



and in the two (gauge invariant) tensor mode polarizations h+ and h x encoded in h*ij (the vector modes disappear 
once the universe becomes isotropic). 

Expression (|B4[) gives the metric perturbations before any gauge is chosen. By equating them with our expressions 
(|Bip wc find how our modes can be decomposed into 3d scalar, vector, and tensor perturbations. We can then use the 
resulting expressions to write 1Z, h+ and h x in terms of our modes. To make this identification, we should spell out 
explicitly how the components of the 3d vector and tensor modes enter in (|B4[) , accounting for their transversality, 
and traceless properties. We do so in the i M coordinate system, for which the momentum of the mode is fc M = 

(O, 0, y/k\ + k\ 2 + 

k-T3 ) ■ I n this system, eqs. (|B4[) give 



5 9uu ( k ) = 



( — 2$* a_B»i ai?*2 aikB* 

a 2 (-2*, + h+) a 2 h y a 2 ikE^ 

a 2 (-2*»-/i + ) a 2 ikE* 2 

\ a 2 (-2** - 2k 2 E„) J 



(B6) 



The entries of this metric can be now identified with those of our metric (|B1[) . transformed according to (|B3I) . We 
obtain 



n 



k 2 T T ff, 
2k 2 + 6 



k 2 

k 2 



kh k^ g 



B r 



(B7) 



Finally, we rewrite our three variables \f r , 5(f), and B v in terms of the canonically normalized modes introduced in 
ID) and (HH (see also eq. ([16]) ). This leads to 



i? 



72 



x 3 / 2 M. 



fl4 



iV2 

j 3 / 2 M T 



(B8) 



As a check, we can verify that the our evolution equations reduce to the standard ones in the limit of isotropic 
background. Using pa = b = in the explicit expressions given in Appendix [A] we find that K s . v — and ft 2 v are 
diagonal in this limit. Therefore, all the canonical modes are decoupled. One can then show that, in this limit the 
equations (l48l) give 



Y±_ 

y/a 



'a 



k z - — 

z 



k 2 - a - 
a 



a 2 
a 







(B9) 



where prime denotes differentiation with respect to conformal time 77, related to the physical time t by drj = dt/a. 
From these expressions, and from eqs. (|B8[) we find 



v"+[k 2 -—\v = , (ah+, x )"+lk 2 J(ofe +1 x)=0 



(BIO) 
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where v = zlZ. These equations are the standard ones [20|, confirming that our formalism reduces to the standard 
one in the limit of isotropic background. 
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